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Abstract 

We present a general study on exclusive semileptonic decays of heavy (B, D, 
B s ) to light (n, p, K, K*) mesons in the framework of effective field theory of 
heavy quark. Transition matrix elements of these decays can be systematically 
characterized by a set of wave functions which are independent of the heavy 
quark mass except for the implicit scale dependence. Form factors for all 
these decays are calculated consistently within the effective theory framework 
using the light cone sum rule method at the leading order of 1/toq expansion. 
The branching ratios of these decays are evaluated, and the heavy and light 
flavor symmetry breaking effects are investigated. We also give comparison of 
our results and the predictions from other approaches, among which are the 
relations proposed recently in the framework of large energy effective theory. 
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I. INTRODUCTION 



Heavy meson decays, both inclusive and exclusive, have long been an interesting subject 
in both experimental and theoretical study. These decays play their special role in extracting 
the Cabibbo-Kabayashi-Maskawa (CKM) matrix elements and probing new physics beyond 
the standard model (SM). The inclusive decays of heavy mesons are more difficult to mea- 
sure but theoretically cleaner than exclusive decays. On the other hand, exclusive decays 
are cleaner in experimental measurements but more difficult in theoretical calculations as 
they require the knowledge of form factors, which contains long distance ingredients and 
has to be estimated via nonperturbative methods such as sum rules, lattice simulations or 
phenomenological models. 

The heavy to light exclusive decays can be grouped as semileptonic decays and rare 
decays. Besides using lattice calculations and quark models [j3[-|T2|l, these decays have 
been analyzed by using sum rules in full QCD [|T3|-pT|]. Since the heavy meson B^ or 
contains one heavy quark and one light quark, it is expected that the heavy quark symmetry 
(HQS) and the effective field theory of heavy quark may help to improve our understanding 



on heavy to light decays. Recently, some work has been done in this direction p2|-p5|. In 
particular, in Refs. p4| , ^5| the B — > -k{p)Iv decays are investigated within the framework 
of heavy quark effective field theory (HQEFT) |2B|-f2S], and the relevant form factors are 
calculated at the leading order of heavy quark expansion, using which \V u b\ is extracted. 

This paper will provide a more general discussion on exclusive heavy to light decays 
within the framework of effective theory. We shall apply the heavy quark expansion (HQE) 
and light cone sum rule (LCSR) techniques to more heavy (B, D, B s ) to light (n, p, K, 
K*) exclusive semileptonic decays. Namely, we extend the study to D decays, and also 
extend it to decays into kaon mesons. We know that the reliability of HQE is different for B 
and D mesons, and that the SU(3) symmetry breaking effects arise in the kaon systems. So, 
through a consistent study on decays of both bottom and charm mesons with the final mesons 
including both nonstrange and strange ones, we aim at a general view on the applicability 
of the combination of HQEFT and LCSR method in studying heavy to light semileptonic 
decays. 

The paper is organized as follows: In section [H], we first formulate the heavy to light 
transition matrix elements in the framework of HQEFT, and then present the light cone 
sum rules for the heavy flavor independent wave functions. Some of the analytic formulae 



are found to be similar to those presented in Refs. [24,25] and so have a general meaning. In 



reviewing them, we pay our main attention to the relations and differences among different 
decay channels. Section [Hj contains our numerical analysis on the heavy to light transition 
form factors. The numerical results are compared with data from other approaches. Based 
on the results of section [Hi] , branching ratios are evaluated and discussed in section |V. 
Finally a short summary is presented in section |V|. 



II. WAVE FUNCTIONS AND LIGHT CONE SUM RULES 

For convenience of discussions, we denote in this paper the light pseudoscalar and vector 
mesons as P and V respectively, and use M to represent the heavy mesons B, B s and D. 
Then the decay matrix elements and form factors can be written in a general form as follows 
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(P(p)\q^Q\M(p + q)) = 2/+(gV + (/+(?*) + 



(2.1) 



<^(P, e*)|g7 M (l - 7 & )Q|M(p + g)) = -i(m M + 2 )e*" + i 



x(2p + qY + i 



A 3 (q 2 



m M + m v 



[e* ■ {p + q))q» + 



2V(q 2 



m M + m v 



m M + m v 



[e*-(p + q)) 



e^e* a {p + q)p Pj , (2.2) 



where q in the currents (not to be confused with the lepton pair momentum) represents light 
quarks (u, d or s), and Q denotes any heavy quark {b or c). mjw, mp^v) are the heavy and 
light pseudoscalar (vector) meson masses, respectively. In HQEFT, the leading order matrix 
elements in the 1/ttiq expansion can be simply expressed as the following trace formulae 
22|-g5| 



(P(p)\qTQ+\M v ) = -Tr[n(v,p)TM v ], 
(V(p,e*)\qTQ+\M v ) = -iTr[Q{v,p)TM v ) 



(2.3) 
(2.4) 



with 



w(v,p) = r/ 5 [A(v ■ p, fi) + fB(v ■ p, //)], 

£l{v,p) = L^v-pjf + L 2 {vp){v ■ e*) + [L 3 (vp)^* + L 4 (v ■ p)(v ■ e* 



(2.5) 
(2.6) 

where p M = p^/v ■ p. in Eqs. fl2.3j ), (|2.4|) is the effective heavy quark field variable 
introduced in the HQEFT pBLB7| , which carries only the residual momentum = p^—mqv^ 
with f M the heavy meson's velocity. Correspondingly M v is the effective heavy meson state. 
It is related to the heavy meson state M in Eqs. ( J2.1|) and (|2.2[ ) by the normalization of the 
hadronic matrix elements ||27|| : 



(ir(p)\qTQ\M) 



1 



A 



-.{(n(p)\qTQt\M v ) + 0(l/m Q )} 



(2.7) 



M 



with Am = rn, m — vnQ. M. v in Eqs.([2~3l 
in the HQEFT (Z7|, 

M, = 



is the heavy pseudoscalar spin wave function 



T l + f 5 



(2.8) 



which exhibits a manifest heavy flavor symmetry. Here A = limM Q _ j . (X) A^ is the heavy flavor 
independent binding energy. 

The relations between the form factors defined in Eqs.( |2.lD , Q2.2D and the universal wave 
functions defined in Eqs.( |2.3| )-( |2.6| ) can be derived straightforwardly. One has 



f±(Q 2 ) 
Ai(q 2 ) 
A 2 (q 2 ) 



m M 



\ -t — i A ■ P) ± B ■ P) } + 

\ A M v-p 



m M A 



rriM + my \ A 
2(m M + m v 



M 



{L x {v p) + L 3 {v ■ p)} + 



m M A i L 2 {v -p) L 3 (v -p) - L A (v-p) ^ _ 

L rv 2 ' o™ I.. „\ J ' 



M 



2m 2 M 



2m M (v-p) 



(2.9) 
(2.10) 
(2.11) 



3 



A 3 (q 2 ) = 2(m M + m v )* 



mMA L 2 (v-p) L 3 (v ■ p) - L 4 (v ■ p) 
\ A M 2m M 2m M (v ■ p) 



m M Am M + m v 
\ A M m M (v-p) 



Note that the form factors introduced in Eqs.( [2.1|) , ( p.2|) are heavy flavor dependent. But 
the functions A, £> and Lj(z = 1,2,3,4) in Eqs. (|2.5|) , (|2.6| ) are leading order wave functions 
in the 1/tuq expansion, so they should be (at least explicitly) independent of the heavy 
quark mass. An advantage using the effective field theory of heavy quark is that it enables 
one to formulate wave functions conveniently in such a heavy quark mass independent way. 
When the wave functions A, B and Lj(i = 1, 2, 3) are estimated to a relatively precise extent, 
then the form factors /+, /_, Ai(i = 1, 2, 3) and V for different decay channels can be easily 
obtained by adopting the relevant parameters such as masses and binding energies of the 
mesons. In other words, Eqs. (|2.9|) - (|2.13|) show the relations among different decays. From 
this point of view, considering the whole group of heavy to light semileptonic decays, one can 
say that to certain order of the I/tuq expansion, the HQS and HQE simplify the theoretical 
analysis and reduce the number of independent functions, though it is well known that for 
an individual decay the number of independent functions does not decrease. 

In light cone sum rule analysis, the wave functions are explored from the study of appro- 
priate two point correlation functions. For example, for decays into pseudoscalar and vector 
light mesons, using the interpolating current Qij 5 q for the pseudoscalar heavy mesons, one 
may consider the functions: 

P^p,q)=i J d A xe iq - x {P(p)\T{q{x)^Q{x),Q{0)rfq{0)}\0), (2.14) 

V»(p,q)=i I ci 4 a;e-^ x (^(p,e*)|T{g(0)7^(l-7 5 )Q(0),Q(a;)z7 5 g(a;)}|0). (2.15) 



In applying for the sum rule method, on the phenomenological considerations, a complete 
set of states with the heavy meson quantum numbers are inserted into the above two point 
functions, i.e., between the two currents. For the insertion of the ground states of heavy 
mesons, one obtains meson pole contributions. While for the insertion of higher resonances, 
the results are generally written in the form of integrals over physical densities pp(v ■ p,s) 
and pv{v-p, s). In the resulting formulae, the matrix elements can be expanded in powers of 
1 / tuq in the effective theory. In this paper we consider only the leading order contributions 
in the heavy quark expansion. Similar to Refs. |23J2H|, we have 



F"(p, q) = 2iF A ^ + BP \ + r ds Pp ( V ' P,S } + subtractions, (2.16) 
2 Am — 2v ■ k J So s — 2v ■ k 

2F e* ■ v 

V»{p, q) = ^ TT-TjU^i + L ^ Hl - W(e* ■ v) - (L 3 - U)p» 

2Am — 2v ■ k v ■ p 

_iAe^ £ W + [°° ds Pv ^ V ' P ' $ + subtractions, (2.17) 
v ■ p J s s — 2v ■ k 

where F is the scaled decay constant of the heavy meson at the leading order of 1/ rriQ |28| . 
Note that when studying B — >• p decay in ||25|| , there is an overall factor multiplied 



to the meson pole contribution. Though such a factor can be obtained in the effective field 
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theory, it can be written in the form 1 + 0(1/ tuq). Since we want Eqs.( |2.16| ) and ([2.17|) 
represent only the leading order contribution in heavy quark expansion, it is consistent to 
miss the factor ™q^ m here. This change may enlarge the B — > plv decay form factors 



obtained in by the same ratio, as can be seen in the next section. 

Heavy quark expansion can be performed to the correlators ( [2.14| ), ( [2.15|) in the frame- 
work of effective theory, after which the field variables and meson states can be replaced by 
their counterparts in effective theory, and those correlators turns into 

P»(p,q) = tJ d 4 xe^- m ^^P(p)\Tq(x)YQt(x),Qt(0)^^ (2-18) 
V»(p,q)=iJ d 4 a;e-^ +m ^ 

+ 0{l/m Q ). (2.19) 

We will include for light pseudoscalar and vector mesons the distribution amplitudes up 
to the same order as in Refs. [24, 25fl. In other words, we consider the ir meson distribution 
amplitudes up to twist 4 and p meson distribution amplitudes up to twist 2. These ir, p 
distribution amplitudes are defined by 

(n( P )\u(x)Yl 5 d(0)\0) = -ijfU [ 1 due^ x [Mu) + x 2 g 1 (u)] 



+ f v (a?--?-) / due iu ^g 2 (u), 
x ■ p Jo 

(7r(p)|n(x)z 7 5 d(0)|0) = f ^ C due iup - x <j>Ju), 

m u + m d Jo 

(n(p)\u(x)a pul 5 d(0)\0)=i(p p x u -p u x p )—^-^- / due*"* 

6{m u + m d )Jo 

(p(p,6*)\u(0)a» u d(x)\0) = -iti(elp v -e: Plx ) [' due^M^), 



e • x 



(p(p,e*)\u(0)^d(x)\0) = f p m pPp / due^^u) 

p ■ x Jo 

+ f P m P (e; - p,—) f due^g^iu), 
p p ■ x Jo 

(p(p,e*)\u{0) W d{x)\0) = \f p m p e pva ^p a x p ^due l ^ x g { l\u). (2.20) 

For K and K* mesons, we need only replace the ir and p mesons and d quark in the lhs. 
of the above equations into K and K* mesons and s quark, and at the same time change 
the quantities and distribution amplitudes related to n and p mesons to those related to K 
and K* mesons in the rhs. of the equations. For decays into Kaon mesons, the light flavor 
SU(3) breaking effects may exhibit via both the light meson related quantities and the Kaon 
meson distribution amplitudes. 

The standard procedure of light cone sum rule method is to calculate the correlation 
functions in deep Euclidean region by using QCD or effective theories, then equate the results 
with the phenomenological representations. In searching for reasonable and stable results, 
the quark-hadron duality and Borel transformation are generally applied to both sides of 
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the equations. Different from the previous LCSR analysis in QCD framework [|13| -|T6|1, in 
the current study, we instead calculate the correlation functions in HQEFT, i.e., we adopt 
the Feynman rules in effective theory. In particular, we use ^ / °° dtS(x — y — vt) for the 
contraction of effective heavy quark fields Q^(x) and Q+(y). As one can see in Refs. fl24U25| , 
the theoretical calculations can often become simpler in the framework of effective theory 
than in QCD. Adopting procedures similar to Refs. p^j25|| , we get 



U(K) 




A P 


Jo 


fn(K) 




4F 


Jo 




d 


6y 




1 r s ° 






dse 


" AF Jo 




= 0, 




1 r s o 






dse 


"AF Jo 




1 r s ° 






dse 


= If Jo 



2y 



y 

[ T y fvmv{ h 



.(«), 



(2.22) 

' { l\u))]u=^ (2.23) 

(2.24) 

fU±(u)] u=ty , (2.25) 



UV) = ^j Q dse^ l -f v m v [^{u) - g ( l\u) + ^#(«))]«=^ (2-26) 

L>2(y) equals zero in the present approximation since no twist 2 distribution amplitudes 
contribute to it. As a consequence, one can see from Eqs. fl2.liD , ( [2.1 2|) that A 2 and A% have 
the same absolute value but opposite signs at the leading order. 



Before proceeding, we would like to address that though the formulae ( |2.21|) - (|2~26|) are 



explicitly independent of the heavy quark mass, it does not mean that the values of these 
functions for decays of different heavy mesons are necessarily the same. This is because that 
one needs to take into account different energy scales for different heavy mesons in calculating 
the corresponding functions. It is this difference that makes the distribution amplitudes and 
other light meson parameters change their values for decays of different heavy mesons. There 
is also an exponent e 2AM ^ T in each formula, which indicates the dependence on the binding 
energy of the heavy meson. Since A# — A^ is small, this exponent only introduces a slight 
difference among the wave functions for B and D decays. On the other hand, there are 
light meson distribution amplitudes, light meson masses and other light meson parameters 
in Eqs.( |2.2l| )-( |2.2tj| ). Consequently, the resulting numerical values for the universal wave 
functions may be different for different light mesons. These heavy and light flavor symmetry 
breaking effects will be explored numerically in the following sections. 



III. NUMERICAL ANALYSIS OF THE FORM FACTORS 

The light cone distribution amplitudes embody the nonperturbative contributions, and 
they are of crutial importance for the precision that light cone sum rules can reach. The 
study of these distribution amplitudes constitutes an important and difficult project. They 
have been studied by several groups. The asymptotic form and the scale dependence of 
these functions are given by perturbative QCD [p9| ,j30[|. 
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For light pseudoscalar mesons, the leading twist distribution amplitude is generally writ- 
ten as an expansion in terms of the Gegenbauer polynomials C^J 2 {x) as follows: 



<t> viK) (u, //) = 6u(l - u)[l + Yl a< K \fi)C^ 2 (2u - 1)]. (3.1) 

n=l 



One should perform the sum rule analysis at appropriate energy scale \x. In the processes of 
Br a ) and D decays, the scales can be set by the typical virtualities of the heavy quarks, for 



example, //& = \Jvnr B — m 2 « 2.4GeV and fx c = yrn? D — m 2 1.3GeV, respectively f31|j . 
For Ti meson, we use HSU 



4(/i c ) = 0.41, <(/i c ) = 0.23, a\ = al = 0, (3.2) 

while for the kaon distribution amplitude 4>k, when the light flavor SU(3) breaking effects 
are considered, we may use fllql 



afOO = 0.17, af (// c ) = 0.21, af (/i c ) = 0.07, af (/x c ) = 0.08, (3.3) 

where the nonvanishing values of the coefficients ai, 03 imply the asymmetric momentum 
distributions for the s and u, d quarks inside the K meson. Since the Gegenbauer moments 
Oj renormalize multiplicatively, the values of di(nb) can be obtained from Eqs. (|3.2|) and (|3.3|) 
through the renormalization group evolution. 

We neglect the SU(3) breaking effects in the twist 3 and 4 distribution amplitudes in- 



cluded in this paper. This is justified by the analyses in Ref. [ 3"2"| , which indicates that these 
breaking effects would influence the light cone sum rules very slightly. Therefore we take for 
both Ti and K the following twist 3 and 4 distribution amplitudes |15|,^I|,^3| : 



<f> p (u) = 1 + -B 2 [3{2u - l) 2 - 1] + -5 4 [35(2u - l) 4 - 30{2u - l) 2 + 3], 
2 8 

(j) a {u) = 6u(l - w){l + lc 2 [5(2u - l) 2 - 1] + ^C 4 [21(2« - l) 4 - U(2u - l) 2 + 1]}, 

2 8 

5 1 6 

gi (u) = -5 2 u 2 {l -uf + -e5 2 [u(l - u){2 + 13u(l - u) + 10w 3 logu(2 - 3m + -u 2 ) 

Z Z 

+ 10(1 - uf log((l -u)(2- 3(1 - u) + jj(l - if) 2 ) )] , 

S 2 ( M ) = y5Vl- M )(2 M -l), (3.4) 
where 

B 2 (/i 6 ) = 0.29, B 4 (// 6 ) = 0.58, B 2 (/i c ) = 0.41, B 4 (// c ) = 0.925, 

C 2 (fib) = 0-059, C 4 (^ 6 ) = 0.034, C 2 (p c ) = 0.087, C 4 (^ c ) = 0.054, 

5 2 (/i b ) = 0.17GeV 2 , 5 2 (/i c ) = 0.19GeV 2 , e(> 6 ) = 0.36, e(// c ) = 0.45. (3.5) 

Besides in these distribution amplitudes, the SU(3) breaking effects also emerge in the 
light meson constants in the coefficients of distribution amplitudes. We take f n = 0.132GeV, 
fx = 0.16GeV and p, n = m 2 /(m u + m d ) with /^(lGeV) = 1.65GeV. For \i K = m 2 K /{m s + 
m u,d)i we use the advocation in Ref. |32| to rely on chiral perturbation theory in the SU(3) 
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limit and so use hk = Z^r- F° r the quantities relevant to heavy hadrons, here we use the data 



evaluated in the previous paper Ref. pSfl . In particular, there we yielded A = 0.53±0.08GeV, 
F = 0.30 ± 0.06GeV 3/2 . 

For decays into light vector mesons, the leading twist distribution functions <p± and </>m 
can also be expanded in Gegenbauer polynomials C^ 2 (x) with the coefficients running with 
the scale and described by the renormalization group method. Explicitly we have 

<f> m {u,n) = 6u(l-u)[l+ £ a^\^)C^(2u-l)}, 

n=2,4,--- 

a^V) = a^")( / , )(^\)(- (ll) -^ (ll) )/(^o) ) (3 . 6) 



where [3 = 11 — (2/3)n/, and r„, rj| are the one loop anomalous dimensions ||3^j35|| . The 
coefficients a^, a\ have been studied extensively in Ref. [0. Here we use the values for p 
and K* mesons presented in that paper, where the SU(3) breaking effects are included for 
K* meson. 

The functions and describe transverse polarizations of quarks in the longitu- 
dinally polarized mesons. They receive contribtuions of both twist 2 and twist 3. In this 
paper we will include only the twist 2 contributions, which are related to the longitudinal 
distribtuion 4>\\(u, p) by Wandzura-Wilczek type relations []35LK7||: 



2 Jo 1 — v Ju v 
g^\u,p) = 2[(l-u) [^v^^ + uf'dv^^]. (3.7) 



o 

The quatities f p and fx* are the decay constants of vector mesons, and and f^* are 
couplings defined via 

(0\ua^q\V(p, e)) = i(e^p u - e v p u )fy (3.8) 
with the light quark q = d or s corresponding to the vector meson V = p or K* . In the 



calculations, we use for these couplings |13|,|36|,p8,39 



f p = 195 ± 7MeV, ff = 160 ± lOMeV, 

f K * = 226 ± 28MeV, = 185 ± lOMeV. (3.9) 



From Eqs.(^-( [2~T3D and (|2~2l|) -(|2~26|), the data for form factors can be obtained by 



using the distribution amplitudes and meson quantities presented above. As an example, 
Fig.l shows the variation of Ai^ p and A^ K * with respect to the Borel parameter T at 
the fixed values of sq = 1.5GeV, 2GeV, 2.5GeV. According to the light cone sum rule 
criterion that both the higher resonance contributions and the contributions from higher 
twist distribution amplitudes should not be too large ( say larger than 30%), our interested 
regions of T are around T = 2GeV for decays, and around T = lGeV for D decays. 
We studied respectively the variations of form factors in all interested decays with respect 
to the Borel parameter T. According to our detailed study, for all form factors, there exist 
reliable regions of Sq and T that well satisfy the requirement of stability in the sum rule 
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analysis. We note that in the original calculation (the first Ref. in [p5fl ) we used a wrong 
sign for the contribution of g± . In the latter study (the second Ref. in [p5|| , or an erratum 
to be published) it is found that the correction of this sign greatly enlarges the value of V, 
and at the same time also improves the stability of the V(T) curves in the sum rule window. 
The threshold energies for all form factors in these decays are found to vary from 0.4GeV 
to 3.5GeV. With the ranges of T and So determined, all the form factors as functions of the 
momentum transfer q 2 can be derived from Eqs. (|2.9| )- (|2.13| ) and (|2.21|) - (|2.26|) . 

It is known that for the final mesons being light pseudoscalars, the light cone expansion 
and the sum rule method will break down at large momentum transfer (numerically as q 2 
approaches near ml) [T5|] . As a result, the curves of wave functions calculated from light 
cone sum rules may become unstable at large q 2 region. Thus at this region we have to 
rely on other approximation such as single or double pole approximation. Here we use for 
B — > 7i transition 

^ - 2 Jr-7k.) (3 - io » 

for large q 2 regions of B decays into light pseudoscalar mesons. And similar monopole 
approximation formulae will be applied to other heavy mesons decaying into light pseu- 
doscalar mesons. In our numerical calculations, we will take: fs* = 0.16 ± 0.03GeV, 



g B - B * = 29 ± 3 @, fago.^ = 2.7 ± 0.8GeV, f D *g D:DK = 3.1 ± 0.6GeV p|, and 
fB*g B *B a K = 3.88 ± 0.31GeV gg. 

As a good approximation, for the behavior of the form factors in the whole kinematically 
accessible region, we use the parametrization 

F(g2) = l-a F q 2 /ml + b F (q 2 /m 2 B ) 2) (3,11) 

where F(q 2 ) can be any of the form factors /+, /_, Ai(i = 1,2,3) and V. Thus each form 
factor will be parametrized by 3 parameters F(0), a F and bp that need to be fitted. For 
the form factor / + , as mentioned above, since the light cone sum rules are most suitable for 
describing the low q 2 region of the form factors and the very high q 2 region is hard to be 
reached by this approach, we shall use the light cone sum rule results at small q 2 region and 
the monopole approximation ( |3.10|) at large q 2 region to fit the 3 parameters in Eq.( |3.11| ). 



For decays into vector mesons, we use only the light cone sum rule predictions in fitting 
these parameters. This is because that the kinematically allowed ranges of B to vector 
meson decays are small compared with the ranges of B to pseudoscalar meson decays, and 
therefore the sum rules are expected to yield reasonable values for most allowed regions of 
q 2 in the vector meson cases. 

At certain values of suitable thresholds Sq, we can fix the parameters for each form factor. 
Our numerical results for those parameters are presented in Table 1-2. For the form factor 
/_, since it is irrelevant to the decay rates when the lepton masses are neglected, we do not 
present their values in these tables. In this paper we calculate them directly from sum rules, 
and the results at large recoil regions will be used in the next section in comparing with the 
form factor relations derived based on the large mass and large energy expansion. 
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F(0) 


dp 


I, 

o F 


s (GeV) 


B - 


-» ttIu 


/+ 


0.35 ±0.06 


1.31 ±0.15 


0.35 ±0.18 


2.3 ±0.6 


B s - 


-> A7z/ 


/+ 


0.47 ±0.09 


1.12 ±0.25 


0.34 ±0.19 


2.7 ±0.8 


D - 




/+ 


0.67 ±0.19 


1.30 ±0.30 


0.68 ±0.38 


0.8 ±0.4 


D - 




/+ 


0.67 ±0.20 


1.65 ±0.43 


1.28 ±0.52 


0.8 ±0.4 



Table 1. Results for the form factor /+ of heavy to light pseudoscalar meson 
decays, sq is the threshold at which the parameters F(0), ap, bp are fitted. 







F(0) 


of 


bp 


so(GeV) 


B -> plu 




0.29 ±0.05 


0.35 ±0.15 


-0.24 ±0.12 


2.1 ±0.6 




0.28 ±0.04 


1.09 ±0.13 


0.20 ±0.18 


^3 


-0.28 ±0.04 


1.09 ±0.13 


0.20 ±0.18 


V 


0.35 ±0.06 


1.32 ±0.12 


0.46 ±0.10 


B s -> K*lv 


A, 


0.28 ±0.05 


0.82 ±0.06 


0.05 ±0.14 


2.5 ±0.5 


A 2 


0.28 ±0.04 


1.48 ±0.05 


0.62 ±0.08 


^3 


-0.28 ±0.04 


1.48 ±0.05 


0.62 ±0.08 


V 


0.35 ±0.05 


1.73 ±0.02 


0.95 ±0.35 


D -> p/z/ 


A, 


0.57 ±0.08 


0.60 ±0.20 


0.51 ± 1.32 


2.0 ±0.5 


A 2 


0.52 ±0.07 


0.66 ±0.20 


-2.03 ± 1.52 


A 3 


-0.52 ±0.07 


0.66 ±0.20 


-2.03 ± 1.52 


V 


0.72 ±0.10 


0.95 ±0.06 


2.60 ±3.62 




Ax 


0.59 ±0.10 


0.58 ±0.12 


0.11 ±0.28 


2.0 ±0.5 


A 2 


0.55 ±0.08 


0.84 ±0.31 


-1.29 ± 1.12 


A 3 


-0.55 ±0.08 


0.84 ±0.31 


-1.29 ± 1.12 


V 


0.80 ±0.10 


0.86 ±0.65 


2.71 ± 1.83 



Table 2. Results for the form factors of heavy to light vector meson decays. 
sq is the threshold at which the parameters F(0), ap, bp are fitted. 



The values of F(0) for B — >• p decay presented in Table 2 are slightly different from those 
given in Refs. p5| . The reason is that the overall factor in Eq.(3.13) of [2l| is missed 
in the current calculation, as mentioned in the previous section. Uncertainties to the form 
factors could arise from the meson constants, the light cone distribution amplitudes and 
the variation of thresholds Sq. We notice that the latter comprises the largest uncertainty, 
which may be 15% or so. Variations of other imput parameters within their allowed ranges 
which we have discussed would increase the uncertainties to about 20%. So, including 
uncertainties from higher twist amplitudes and other systematic uncertainties in light cone 
sum rule method, we quote an uncertainty of about 25%. 

Besides the SU(3) symmetry breaking effects arise from relevant light meson parameters, 
it is found in our investigation that SU(3) symmetry breaking effects considered in the 
Gegenbauer polynomial moments a,i(i = 1,2,3,4) invoke the changes of F(0) by no more 
than 15%. 

The form factors as functions of q 2 are shown in Fig. 2. Note that in the figures some 
curves almost coincide with each other, which indicates that the SU(3) breaking effects are 
small in those cases. 
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Table 3 is a comparison of the form factor values at q 2 = obtained in this work and 
those obtained from other approaches, including quark models, lattice calculations, and also 
sum rules in QCD framework. 



B -> n(p)lu 



Ref. 

this work 
QM § 
QM § 
Lat || 
Lat H 

SR 



17, IS 



SR p| 



/+(0) 
0.35 ±0.06 

0.29 

0.27 
0.50(14)1£ 
0.28 ±0.04 

0.305 
0.28 ±0.05 



^i(O) 
0.29 ±0.05 

0.26 

0.26 
0.16(4)1 



\+22 
-16 



0.26 ±0.04 



A 2 (0) 
0.28 ±0.04 

0.24 

0.24 
0.72(35)1^° 



0.22 ±0.03 



0.35 ±0.06 

0.31 

0.35 
0.61(23)1;? 

0.34 ±0.05 



B s -> K(K*)lu 



this work 
QM § 



0.47 ±0.09 
0.31 



0.28 ±0.05 
0.29 



0.28 ±0.04 
0.26 



0.35 ±0.05 
0.38 



D -> n(p)lu 



this work 
QM § 

qm i 

QM 

Lat || 
Lat § 
Lat M 



SR p| 



SR |21| 



0.67 ±0.19 

0.69 

0.67 

0.69 
0.68(13)1^° 
0.64(5)18? 
0.65 ±0.10 
0.65 ±0.11 
0.50 ±0.15 



0.57 ±0.08 
0.59 
0.58 
0.78 
0.59(7)i| 

0.65 ±0.07 
0.5 ±0.2 



0.52 ±0.07 

0.49 

0.42 

0.92 
0.83(20)1^ 2 

0.55 ±0.10 

0.4 ±0.2 



0.72 ±0.10 

0.90 

0.93 

1.23 
1.31(25)111 

1.1 ±0.2 

1.0 ±0.2 



D -> K(K*)lu 



this work 



E791 g| 
QM § 

qm i 

Lat § 
Lat || 
SR P| 



SR |T| 



0.67 ±0.20 

0.78 

0.78 
0.71(12)1^° 
0.73 ±0.07 
0.76 ±0.03 
0.60 ±0.15 



0.59 ±0.10 
0.58 ±0.03 

0.66 

0.66 
0.61(6)1? 
0.70 ±0.07 

0.50 ±0.15 



0.55 ±0.08 
0.41 ±0.06 

0.49 

0.43 
0.83(20)1^ 2 
0.6 ±0.1 



0.60 ±0.15 



0.80 ±0.10 
1.06 ±0.09 

1.03 

1.04 
1.34(24)112 
1.2 ±0.2 

1.10 ±0.25 



Table 3. Comparison of the results for semileptonic decay form factors (at 
q 2 = 0) from different approaches (QM: quark model; Lat: lattice calculation; 
SR: sum rules in QCD framework; E791: data extracted from experimental 
measurements) . 



The form factor ratios R v = V(0)/Ai(0), R 2 = A 2 (0)/A 1 (0) for D -> K*lv decay have 
been measured by several groups. We present in Table 4 the comparison of our results 
for these ratios with the measurements. Our result for Ry agrees well with the latest 
measurement HI] but smaller than other measurements. R 2 obtained in this work is also in 



good agreement with the BEATRICE measurement [41 
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R 2 


A 1 (0) 


this work 


1.36 ±0.39 


0.93 ±0.29 


0.59 ±0.10 


BEA1R1UE |41 


1.45 ± 0.23 ± 0.07 


1 nn i n 1 r i n no 

1.00 ± 0.15 ± 0.03 




E791 H 


1.84 ±0.11 ±0.08 


0.71 ±0.08 ±0.09 


0.58 ±0.03 


E687 || 


1.74 ±0.27 ±0.28 


0.78 ±0.18 ±0.10 




E653 g 


2.00^ ±0.16 


± 0.11 




E691 g§ 


2.0 ±0.6 ±0.3 


0.0 ±0.5 ±0.2 





Table 4. Comparison of measurements and theoretical predictions for form 
factor ratios and Ai(0). 



With the consideration that in heavy to light decays the final light meson usually carries a 
large recoil momentum, the large energy effective theory (LEET) [^ -501 has been proposed 
to study heavy to light transitions in the region of large energy release. In that framework 
and neglecting QCD short distance corrections, it is quoted that there are only three form 
factors needed describing all the pseudoscalar to pseudoscalar or vector transition matrix 
elements in the leading order of heavy quark mass and large energy expansion. In particular, 
taking into acount contributions of second order in the ratio of the light meson mass to the 
large recoil energy, the recent work [|5(J derived interesting relations concerning form factors 
for semileptonic B decays to light pseudoscalar and vector mesons: 



U(q 2 ) 



m M 
2E F 



1 + 



ml 



m 2 M 



)fo(q 2 ), 



(3.12) 



m M 



El 



mf 



m M + my 



rriM + my 
2E~ F 

m M 



(1 _ 2^ 

m M + m v m 2 M 
where Ep is the on-shell energy of the final light meson: 



Hip 



(3.13) 



E, 



m M + m 



P(V) 



Q 



2m,M 

fo(q 2 ) and A (q 2 ) are form factors defined in the way of Ref. 
represented by the form factors defined in this paper as follows, 

fo(q 2 )= / 2 f-{q 2 ) + U{q 2 ), 



(3.14) 

and they can be 



m 2 M — m 2 p 



A (q 2 



1 



2mi 



-[(m M + m v )A 1 (q 2 ) - (m M - m v )A 2 (q 2 



m M + my 



Mq% 



As a result, relations (|3.12|) -( [3.13|) are equivalent to 

;fo(q 2 ) 



U(q 2 ) 

V(q 2 ) 
A 2 (q 2 ) 



m 2 M + m 2 P 



m 2 M ± m 2 P - q 2 
(m M + m v ) 2 



2m 



M 



E F 



-A(q 2 



m 



v 



m 2 M q 2 



2m v q 2 



^2-^A 3 (q 2 

™>M<r 



(3.15) 
(3.16) 

(3.17) 
(3.18) 

(3.19) 
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Note that in the heavy quark and large recoil expansion, the leading order contribution of 
the denominator in Eq. (|3.19| ) is just —m 2 M q 2 A 2 (q 2 ). As a result, one gets A 2 (q 2 ) = —A 2 (q 2 ) 
at the leading order of LEET, which is consistent with our calculation in this work (see 
section [D]) . 

One can now have a numerical comparison between our results for form factors with the 
LEET predictions. In Figs. 3-5, we show the form factor ratios f+(q 2 )/ fo(q 2 ), V(q 2 )/Ai(q 2 ) 
and A 2 (q 2 ) j A 3 (q 2 ) derived from our sum rule calculations and from the relations ( |3.17| )- 
03.19| ). In Fig. 3, the ratio f+{q 2 )/ fo{q 2 ) from our calculation and that from the LEET 
relation ( |3.17| ) are close to each other at regions near the maximum recoil point, but they 
quickly split from each other as q 2 increases. Similar variations can be observed for the B^ 
decay ratio V(q 2 )/Ai(q 2 ) in Fig.4(a)-(b). For D decays to light vector mesons, however, 
in Fig.4(c)-(d) large difference exists between our results and the predictions of currently 
known LEET relation (|3.18|) . In Fig.5(a)-(b) (i.e., for decays), the ratio A 2 (q 2 ) / A 3 (q 2 ) 
of our results and the LEET prediction are compatible in almost the whole kinematically 
allowed regions except the region very near to q 2 = point. In Fig.5(c)-(d), our results for 
the D decay ratio A 2 (q 2 ) / A 3 (q 3 ) and the LEET relation ( |3TT9D are again incompatible. 

As a general view, our results are compatible with the relations (|3.17| )- (|3.19|) at appro- 
priate regions for B^ decays. On the other hand, to the expansion order considered here, 
discrepancy occurs between the two methods in studying D decays. One reason for this may 
be of course that the charm quark (hadron) mass is not very large, and therefore large 1 / m c 
corrections both to our results of leading order in HQE and to the LEET approximations 
(p.l7p-(pTT^) are possible. This should be clarified by the future study. 



IV. BRANCHING RATIOS 

With the form factors derived in the previous sections, one can extract the values of 
relevant CKM matrix elements from the experimental measurements of branching ratios. In 
Refs. [24 1 and [25], \V u b\ is extracted in this way from B — > 7r(p)/z/ decays with using the 



branching ratio measurements. The results obtained there were \V u b\ = (3.4±0.5±0.5) x 10 -3 
and \ V u b\ = (3.7 ± 0.6 ± 0.7) x 10~ 3 via the two decays, respectively, where the first (second) 
errors correspond to the experimental (theoretical) uncertainties. 

On the other hand, the decay widths and branching ratios of exclusive semileptonic 
decays can be predicted if we know the values of relevant CKM matrix elements. When the 
lepton masses are neglected, we have for decays to pseudoscalar mesons: 

dT _ G 2 F \V qQ \ 2 m 2 M + m 2 P - q 2 2 2\3/2 {f ( 2^2 uu 
dq 2 ~ 24vr3 U 2m M } ™ pj [U{Q jJ ' ^ 

and for decays to vector mesons: 

j-p r*2 it/ 12 

T-2 = ^r^r^ViHo + H l + H-) ( 4 - 2 ) 

dq 2 19271*171%.; 

with A = ijn 2 M + my — q 2 ) 2 — 4m 2 M my and the helicity amplitudes defined as 

A 1 / 2 



H± = (m M + m v )A 1 {q 2 ) T ; V(q 2 ) 
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H ° = o 7t{(^m -4- <? 2 )(™a/ + myj^fe 2 ) A 2 (g 2 )}. (4.3) 

In this paper we would like to predict the branching ratios from the values of CKM 
matrix elements: \V u b\ = 0.0037, \V c d\ = 0.22, \V CS \ = 0.97. Finishing the integration over 
q 2 , we obtain the widths and branching ratios in Table 5. In that table we also list the 
experimental measurements of branching ratios given by Ref. 





r (in I^qI^s-^ 


Br 


Br ( measurement s) 


B -> ttIv 


10.2 ± 1.5 


(2.1 ±0.5) x 10" 4 


(1.8 ±0.6) x 10" 4 


B s -> Klu 


14.3 ±2.8 


(2.9 ±0.7) x 10- 4 




B -> p/z/ 


15.2 ±4.2 


(3.2 ± 1.0) x 10~ 4 


(2.6i^;?) x 10" 4 




20.2 ±4.3 


(4.1 ± 1.0) x 10" 4 




D -> tt/z/ 


0.152 ±0.042 


(3.0 ±0.9) x 10~ 3 


(3.7 ±0.6) x 10- 3 


£> -> lflz/ 


0.101 ±0.030 


(3.9 ± 1.2) x 10~ 2 


(3.47 ±0.17) x 10" 2 


D -> plu 


0.071 ±0.015 


(1.4 ±0.3) x 10" 3 




D -> If* Zi/ 


0.050 ±0.014 


(2.0 ±0.5) x 10~ 2 


(2.02 ±0.33) x 10" 2 



Table 5. Decay widths F (in unit |V^q| 2 ps _1 ) and branching ratios (Br) for 
heavy to light meson decays. In deriving the branching ratios we used | V^f, | = 
0.0037, \Vcd\ = 0.22, \V CS \ = 0.97 and the lifetimes of heavy mesons: r^o = 
1.56 ± 0.06ps, r D o = 0.4126 ± 0.0028ps, t Bs = 1.493 ± 0.062ps. 



The largest uncertainty for the branching ratios in Table 5 is about 30%. So, the precision 
of heavy meson decay measurements expected in the new B factories can not be matched well 
unless these theoretical uncertainties could be reduced. This reduction may be reached by 
consideration of both higher twist contributions, better determination of the meson constants 
and the higher order contributions in the heavy quark expansion, which should be discussed 
in the future work. 

Table 6 is a comparison of the results for the ratios Tl/Tt, r + /T_ and the total decay 
rates T for heavy to light vector decays, where T L (T T ) and r + (/T_) represent partial widths 
for longitudinal (transverse) polarization and positive (negative) helicity, respectively. These 
widths and ratios have been predicted via other approaches. In particular, the ratios for 
D — > K*lu decay are available from recent experiments. 
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r L /r 



T 



r + /r_ 



Ref. 



B -> p/i/ 



0.85 ±0.09 

0.82 
0.88 ±0.08 

o.so±85S 

0.52 ±0.08 
0.06 ±0.02 



0.04 ±0.02 



0.007 ±0.004 



15.2 ±4.2 

19.1 
15.8 ±2.3 
13 ± 12 
16-5±1;| 
13.5 ±4.0 
12±4 



this work 
QM @ 

QM |T 
Lat § 
Lat 
SR ET 



SR 21 



K*lv 



0.79 ±0.08 



0.07 ±0.02 



20.2 ±4.3 



this work 



D -> pZi/ 



1.17 ±0.09 
1.16 
0.67 



1.31 ±0.11 



0.29 ±0.13 



0.24 ±0.03 



0.071 ±0.015 
0.087 
0.025 
12.17 ±4.06 
10.15 ±4.67 
0.024 ±0.007 



this work 
QM § 

QM [0 

Lat 
Lat @ 
SR |2l 



1.15 ±0.10 
1.09 ±0.10 ±0.02 
1.20 ±0.13 ±0.13 
1.18 ±0.18 ±0.08 

i.s±8:i ± o.3 
o.6 ± o.3±S:? 

1.28 

1.33 
1.2 ±0.3 
1.1 ±0.2 
0.86 ±0.06 



0.32 ±0.13 
0.28 ±0.05 ±0.02 



0.050 ±0.014 



this work 
BEATRICE EI 



E687 @ 



0.063 
0.058 
0.073 ±0.019 
6.0l?| 



E653 H 

E691 g§] 

WA82 [0 
QM § 

QM i 

Lat 
Lat M 



3.8 ± 1.5 



SR 21 



Table 6. Theoretical predictions and measurements of the ratios Tl/Tt, 
r + /r_ and the decay rates T (in unit |V^q| 2 ps _1 ). The data from BEAT- 
RICE, E687, E653, E691 and WA82 are experimental measurements. 

Similarly, Table 7 presents a comparison of results for heavy to light pseudoscalar decay 
rates. 



Ref. 


B -> nlu 


B s -> Klv 


D -> nlu 


D -> Klu 


this work 

QM I 
QM i 

Lat [0 
Lat [| 
Lat [0 
SR (T9fl 
SR 21| 


10.2 ± 1.5 

10.0 

8±4 
12 ±8 

8.5l?;| 

7.3 ±2.5 
5.1 ± 1.1 


14.3 ±2.8 


0.152 ±0.042 

0.196 

0.165 
0.162 ±0.041 
0.114 ±0.073 

0.13 ±0.05 
0.080 ±0.017 


0.101 ±0.030 

0.102 

0.101 
0.096 ±0.021 
0.072 ± 0.036 

0.094 ±0.036 
0.068 ±0.014 



Table 7. Rates in unit |lgQ| 2 ps 1 for semileptonic decays to pseudoscalar 
mesons from different approaches. 
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V. SUMMARY 



In summary, we have studied the exclusive semileptonic decays of heavy to light mesons. 
The form factors for the decays B(D, B s ) — ► n(p, K, K*)lv have been consistently calculated 
by using the light cone sum rule method in the effective field theory of heavy quark. The 
HQS leads to simplification in studying heavy to light tansitions as to certain order of the 
1 / mq expansion the decays of different heavy hadrons such as B and D can be characterized 
by the same set of wave functions which are explicitly independent of the heavy quark mass, 
though the HQS does not reduce the number of independent form factors needed for an 
individual decay. In such calculations, the uncertainties for the form factors are generally 
about 25%, which, together with the meson constants, may give the branching ratios with 
a total uncertainty up to 30%. We have also estimated the light flavor SU(3) symmetry 
breaking effects in these semileptonic decays and found that those effects may influence 
the form factors up to a total amount of 15%. Our results have been compared with data 
from experiments and from other theoretical approaches. In particular, we have checked the 
compatibility between our results and the form factor relations derived in Ref. []5(| using the 
heavy quark and large recoil expansion. We conclude that the form factors and branching 
ratios of those heavy to light meson semileptonic decays can be calculated consistently based 
on the light cone sum rule approach within the framework of HQEFT. Nevertheless, in order 
to match the expected more precise experimental measurements at B factories in the near 
future, more accurate calculations for the exclusive heavy to light meson semileptonic decays 
are needed. The uncertainties in our results and the discrepancy existing between our results 
and the relations (|3.17| )- (|3.19| ) indicate that the 1/rriQ corrections may be important in charm 



meson decays. This should be investigated in the future work. 

In this paper we have only considered the semileptonic heavy to light decays, but what we 
used here is a general approach, and similar procedures could be applied straightforwardly 
to heavy to light rare decays such as B — > Kll, B — > K*ll and B — > K*j. For such rare 
decays, this framework becomes more interesting because it reduces the number of form 
factors for an individual decay. We will discuss that in another paper. 
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FIGURES 




2.5 3 



(a) (b) 

Fig.l. The form factors A^ p ((a)) and Af^ K * ((b)) as functions of 
the Borel parameter T for different values of the continuum threshold 
so- The dashed, solid and dotted curves correspond to sq =1.5, 2 and 
2.5 GeV respectively. Considered here is at the momentum transfer 
g 2 = QGeV 2 . 
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Fig. 2. Results for the heavy-to-light decay form factors from light cone 
sum rule study. The solid, dashed, dot-dashed and dotted curves corre- 
spond to B -► 7r(p), £ s -► X (if*), Z? -► 7r(p) and D if (if*) decays, 
respectively. The dot-dashed and dotted curves in figure (a) almost co- 
incide with each other, and so do the solid and dashed curves in figures 
(b)-(d). As(q 2 ) is not shown here as one has A%{q 2 ) = —A2(q 2 ) at the 
leading order considered. 
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Fig. 3. Results for the ratio f+{q 2 )/ fo(q 2 )- (a), (b), (c), (d) are for 
B — ► irlv, B s — ► iTZi/, D — > -zr/^, D — > ET/^ respectively. Solid lines 
are our sum rul e res ults, while the dashed lines are produced from the 
LEET relation (|3~17|) . 
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Fig.4. Results for the ratio V (q 2 ) / Ai(q 2 ) . (a), (b), (c), (d) are for 
B — ► /)/^, i? s — ► K*lv, D — ► p/^, D — > respectively. Solid lines 

are our sum rul e res ults, while the dashed lines are produced from the 
LEET relation (|3.18[) . 
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Fig.5. Results for the ratio A 2 (q 2 )/A 3 (q 2 ). (a), (b), (c), (d) are for 
B — > plv, B s — * K*lv, D — ► plu, D — > respectively. Solid lines 

are our sum rul e res ults, while the dashed lines are produced from the 
LEET relation (|3.19[) . 



23 



